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Abstract 



Employing the Noether charge technique and Visser's Euclidean approach the en- 
tropy of the nonlinear black hole described by the perturbative solution of the system of 
coupled equations of the quadratic gravity and nonlinear electrodynamics is constructed. 
The solution is parametrized by the exact location of the event horizon and charge. Spe- 
cial emphasis in put on the extremal configuration. Consequences of the second choice 
5^ ' of the boundary conditions, in which the solution is paramerized by the charge and the 

total mass as seen by a distant observer is briefly examined. 

1 Introduction 

m 

Recently, a great deal of efforts have been devoted to the important issue of regular black 
holes. One of the most intriguing solutions of this type have been constructed by Ayon- 
CnJ ' Beato and Garcia [1] and by Bronnikov [2]. In both cases, the line element is a solution of 

the coupled system of equations of nonlinear electrodynamics and gravity. (We shall refer 
to the solutions of this type as ABGB geometries). The former solution describes a regular, 
static and spherically symmetric configuration with the electric charge, Q e , whereas the latter 
one describes a similar geometry characterized by the mass and the magnetic charge Q. For 
certain values of the parameters both solutions describe black holes. On the other hand, 
the no-go theorem proved in Ref. [3] (see also [2, 4]) forbids, for the class of electromagnetic 
Lagrangians with a Maxwell asymptotic in a weak field limit, existence of the electrically 
charged, static and spherically- symmetric solutions with the regular center. It should be 
noted, however, that the electric solution is not in conflict with the non existence theorem, 
as the formulation of the nonlinear electrodynamics [5] employed by Ayon-Beato and Garcia 
(P framework in the nomenclature of Ref. [2] ) differs from the one to which one refers in the 
assumptions of the no-go theorem. Indeed, the solution of Ayon-Beato and Garcia has been 
constructed in a formulation of the nonlinear electrodynamics obtained from the original one 
(F framework) by means of a Legendre transformation (see Ref. [2] for details). Moreover, 
the no-go theorem does not forbid existence of the solutions with magnetic charge as well as 
some hybrid configurations in which the electric field does not extend to the central region. 

The status of the nonlinear electrodynamics in the model considered here is to provide 
a static matter source, perhaps the exotic one, to the field equations. That means that 
the casual structure of the spacetime is still governed by the null geodesies or "ordinary" 
photons rather than the photons of the nonlinear theory. Actually, the latter move along 
the geodesies of the effective space [6, 7]. Outside the event horizon the solution of the 
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ABGB-type closely resembles the Reissner-Nordstrom (RN) geometry both in its global and 
local structure. Important differences appear near the extremality limit. Consequently, the 
Penrose diagrams of the ABGB solution are similar to those constructed for the Rissner- 
Nordstrom solution, with the one notable distinction: instead of the singularity at r = now 
we have the regular interior. 

An attractive feature of the ABGB solutions is possibility to express the location of the 
horizons in terms of the Lambert special functions [8, 9]. Similarly, the Lambert functions [10, 
11] may be used in the discussion of the extremal configurations [12]. 

According to our present understanding a proper description of the gravitational phenom- 
ena should be given by the quantum gravity, being perhaps a part of a more fundamental 
theory. And although at the present stage we have no clear idea how this theory looks 
like, we expect that the action functional describing its low-energy approximation should 
consist of the higher order terms constructed from the curvature tensor, its contractions 
and covariant derivatives to some required order. Among various generalizations of the 
Einstein- Hilbcrt action a special role is played by the quadratic gravity (see for example 
Refs. [13, 14, 15, 16, 17, 18, 19, 20, 21]). Motivations for introducing such terms into the 
action functional are numerous. When invented, for example, the equations of quadratic 
gravity have been treated as an exact formulation of the theory of gravitation. On the other 
hand, it may be considered, quite naturally, as truncation of series expansion of the action 
of the more general theory. Such terms appear generically in the one-loop calculations of the 
quantum field theory in curved background [22]. Moreover, from the point of view of the 
semi-classical gravity, the quadratic terms in the field equations might be treated as some sort 
of the simplified stress-energy tensor. Such a toy model of the renormalized stress- energy 
tensor allows to mimic the fairly more complex sources in a relatively simple way. This ap- 
proach is especially useful when the general pattern that lies behind the calculations of both 
types is essentially the same. Thus, some general features of the full semi- classical solutions 
can be analyzed and understand without referring to otherwise intractable equations. 

It should be noted that any higher curvature theory contain solutions which are unavail- 
able to the theory based on the classical Einstein- Hilbert Lagrangian. This can most easily 
be seen by counting the degrees of freedom: the quadratic gravity is known to posses 8 de- 
grees of freedom whereas the General Relativity has only 2. Moreover, there are solutions 
that are not analytic in the coupling constants, i. e., they do not reduce to solutions of the 
classical Einstein field equations. (For a comprehensive discussion see for example [23] and 
the references cited therein) . Unfortunately, because of complexity of the equations of the 
quadratic gravity it is practically impossible to construct their exact solutions and one is 
forced to refer either to approximations or to numerical methods. The natural method to 
obtain reasonable results consists of treating the higher curvature contributions perturba- 
tively. This approach also guarantees that the black hole exists as the perturbative solution 
of the higher-order solution provided it exists classically [24]. Finally, observe that in the 
perturbative approach the casual structure is determined by the classical metric, however, 
the equations of motion of test particles and various characteristics of the solution acquire 
the first order correction. 

Analyses of the spherically-symmetric and static solutions to the higher derivative theory 
has been carried out in [14, 25, 26, 27, 28, 29, 30]. Specifically, in Ref. [30] the perturbative 
solutions of the ABGB-type to the equations of the effective quadratic gravity have been 
constructed and discussed. In this paper we shall calculate the entropy of such black holes 
using Wald's approach [31, 32, 33] and confirm the final results employing computationally 
independent but closely related Euclidean techniques propounded by Visser [34, 35, 36]. 

2 Basic equations 
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The coupled system of the nonlinear electrodynamics and the quadratic gravity considered 
in this paper is described by the (Lorentzian) action 

S= lL j ' ( R + aR2+ P R abR ah + iRabcdR^ - C (F)) d 4 X, (1) 

where C (F) is some functional of F = F a t,F ab (its exact form will be given later) and 
all symbols have their usual meaning. The cosmological constant is assumed to be zero. 
To simplify our discussion from the very beginning we shall relegate the term involving the 
Kretschmann scalar, R a bcdR abcd , from the total action employing the Gauss-Bonnet invariant. 
The coupling constants a and (3 have the dimension of length squared and throughout the 
paper we shall assume 

ot 8 

where L is the local curvature scale. Assumption that the mass scales associated with the 
linearized equations are real may place additional constrains [20, 37, 38] on a and p. Here, 
however, we shall treat them as small and of comparable order but arbitrary. 

The entropy of the black hole may be calculated using various methods. It seems, however, 
that Wald's technique is especially well suited for calculations in the higher curvature theories. 
Here we shall follow this very approach. Other competing techniques are the method based 
on the field redefinition [33, 39] and Visser's Euclidean approach. 

For the Lagrangian involving the Riemann tensor and its symmetric derivatives up some 
finite order, say n, Wald's Noether charge entropy may be compactly written in the form 
[31, 32, 33] 

S = -2tt [ d 2 x (h) 1/2 J2 (-l) m V {ei .V em) Z^-^ abcd e ab e cd , (3) 

where 

rve\...e m \abcd UUj 



~~ 9V (ei ...V em) i? a6cd ' (4) 

h is the determinant of the induced metric, e a b is the binormal to the bifurcation sphere, and 
the integration is carried out across the bifurcation surface. Actually S can be evaluated not 
only on the bifurcation surface but on an arbitrary cross-section of the Killing horizon. Since 
tabbed = 9ad9bc — (JacQbdj where g ac is the metric in the subspace normal to cross section on 
which the entropy is calculated, one can rewrite Eq. (3) in the form 

/n 
d 2 x h 1 ' 2 £ V (ei ... V em) Z e ^ abcd g ac g bd . (5) 

m=0 

The tensor g a b is related to V a = K a /\\K\\ (K a is the timelike Killing vector) and the unit 
normal n a by the formula g a b = V a Vb + n a rib. 

The general expression describing entropy (5) has been applied in numerous cases, mostly 
for the Lagrangians that are independent of covariant derivatives of the Riemann tensor and 
its contractions. In Ref. [40], however, Eq. (5) has been employed in calculations of the 
entropy of the quantum-corrected black hole when the source term is described by the stress- 
energy tensor of the quantized fields in a large mass limit. Such a tensor is purely geometrical 
and besides ordinary higher curvature terms it involves also RV a V a R and i? a fcV c V c i? ob . 

On the other hand, one can follow an approach propounded by Visser [34, 35, 36]. The 
general formula for the entropy of the stationary black hole with the Hawking temperature 
Tjj is given by 

S =J + ^ Jj P L-L E )K a dZ a + f^V a d^a, (6) 
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where A is the area of the event horizon, s is the entropy density associated with the fluctu- 
ations (ignored in this paper) and finally pl and Le are, respectively, the Lorentzian energy 
density and the Euclideanized Lagrangian of the matter fields surrounding the black hole. 
(All higher curvature terms have been inserted into the Lagrangian describing matter fields.) 
For the specific case of the Einstein-Hilbert action augmented with the higher curvature terms 
(but not covariant derivatives of curvature) Visser's result is equivalent to Wald's formula. 

The coupled system of differential equations describing nonlinear electrodynamics in 
quadratic gravity can be obtained from the variational principle. 

Simple calculations indicate that the tensor F ab and its dual *F ab 1 satisfy the equations 

V.(^f-)-H (7) 

V a *F ab = 0, (8) 

respectively. Differentiating functionally the total action S with respect to the metric tensor 
one obtains equations of the quadratic gravity in the form 

L ab = Qab _ a jab _ pjab = g^ab^ (g) 

jab = 2V h V a J R - 2RR ab + ^g ab (R 2 - 4V c V c i?) , (10) 

jab = v b v a R _ VcV c R ab _ 2 R cd R cbda + l.g ab (R cd R cd - V C V°R) (11) 

nb _ 1 / dC (F) . 1 , 



where 



and 



T a = ^[^^F Cb --/ a C(F)). (12) 



In this paper we shall concentrate on the static and spherically-symmetric configurations 
described by the line element of the form 

rlr 2 

ds 2 = _ e 2Hr) f{r)dt 2 + ^- )+ r 2 dn 2 , (13) 

where 

m-i-^. (») 

r 

The spherical symmetry places restrictions on the components of F a b tensor and, conse- 
quently, its only nonvanishing components compatible with the assumed symmetry are F i 
and F23. Simple calculations yield 

F 23 = Qsin0 (15) 

and 

r*e-^^P-F w =Q e , (16) 

where Q and Q e are the integration constants interpreted as the magnetic and electric charge, 
respectively. 

Since the no-go theorem forbids existence of the regular solutions with Q e ^ in the 
latter we shall assume that the electric charge vanishes. Now, since F = 2F23F 23 , one has 

The stress-energy tensor (12) calculated for this configuration is 

t* = t; = ~c{F) (is) 
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and 

1 dC(F) Q 2 _ _L 
4^ dF ~ ~ 16tt 

Further considerations require specification of the Lagrangian C (F) . Following Ayon 
Bcato, Garcia and Bronnikov let us chose it in the form 



C(F) = F 

where 




(20) 



and & is a free parameter. Inserting Eq. (17) into (20) and making use of Eq. (21) one obtains 

2Q 2 A ,^ 2 Q 2 



r 



£ (^) = -3-( 1 -tenh^). (22) 



The system of coupled differential equations of the quadratic gravity with the source 
term given by (18) and (19) with (22) is rather complicated and cannot be solved exactly. 
Fortunately, since the coupling constants a and are expected to be small in a sense of 
Eq. 2, one can treat the system of the differential equations perturbatively, with the classical 
solution of the Einstein field equation taken as the zeroth- order approximation. Successive 
perturbations are therefore solutions of the chain of the differential equations of ascending 
complexity [41, 42, 43, 44]. It should be noted, however, that the higher order equations are 
probably intractable analytically and the technical difficulties may limit the calculations to 
the first order. 

In the next section, we shall employ perturbative techniques to construct the approximate 
solution to the equations of the quadratic gravity with the source term being the stress-energy 
tensor of the Bronnikov type. Such an approach is expected to yield reasonable results and 
because of complexity of the differential equations, it may be the only way to deal with this 
problem. 



3 Solutions 

To keep control of the order of terms in complicated series expansions we shall introduce a 
dimcnsionless parameter s substituting a^ea and (3 — > e(3. We shall put e = 1 at the final 
stage of calculations. Of functions M (r) and if) (r) we assume that they can be expanded in 
powers of the auxiliary parameter as 

M (r) = M (r) + eM 1 (r) + O (e 2 ) (23) 

and 

V> (r) = eVi (r) + O (e 2 ) . (24) 

First, consider the left hand side of Eq. (9) calculated for the line element (13) with 
the functions M(r) and tp(r) given by (23) and (24), respectively. Making use of the above 
expansions and subsequently collecting the terms with the like powers of e, after some rear- 
rangements, one obtains [30] 

L\ = ~{M' + eM[-eSt), (25) 
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where 



t _ ( 2 M' 8M M; 2M' 2 2M; 5 M Mq M^Mg 
&t — P — o 5 1 5 1 5 

+ + M (3) - _ Mq M< 3 > + r M< 4 > - 2 M M< 4 ^ 

_ q ( 24M M{, _ _ AM? + 8M£ _ 18M Aff _ M „ 2 



+ 2M o^ - 4 M< 3 > + 6M ° r M " 3) + 2 Mq Jlf <»> - 2 r + 4 M M^ ) 

(26) 

and Mq, Mq and Mq^ for i > 3 denote first, second and i— th derivatives with respect to the 
radial coordinate. On the other hand, a simple combination of the components of L h a tensor 

L r r -L\ = (27) 

can be easily integrated to yield [30] 

Mr) = (2« + /3)M (3) - -i(3a + /3)Mq + Ci, (28) 

where C\ is the integration constant. It should be noted that contrary to the case of coupled 
system of the Maxwell equations and quadratic gravity considered in Refs. [25, 26, 27, 29], 
now we have explicit dependence on the parameter a. A comment is in order here regarding 
the independence of the final result calculated for the Maxwell source on the parameter a. 
First, observe that the stress-energy tensor of the electromagnetic field for the spherically- 
symmetric an static configuration with a total charge e assumes simple form 

^ = _ 8^ diag[1 ' 1 '" 1 '" 1] ' (29) 
Therefore, the zeroth-order solution to the (0,0) -component of the equation (9) can be 
written in the form 

M (r) = ~+C, (30) 
2r 

where C is the integration constant. Now, substituting (30) into (26) and (28) it can easily 
be demonstrated that the expression in the second bracket in its right hand side of Eq. (26) 
as well as the expression Mq — QM^/r 2 in (28) vanish. 

One expects that all characteristics of the black hole, such as the location of the horizons 
and temperature could also be calculated perturbativcly. In the latter, for simplicity, we shall 
refer to the perturbative solutions of the quadratic gravity using the names of their classical 
counterparts (the zeroth-order solutions) whenever it will not lead to confusion. 

To develop the model further one has to determine the integration constants and the free 
parameter b. There are, in general, two interesting and physically motivated choices. One 
can relate the integration constant with the exact location of the event horizon, r + , and this 
can easily be done with the aid of the equation 

M(r + )= r -±. (31) 

On the other hand it is possible to express solutions of the system of differential equations 
consisting of (0, 0) component of Eqs. ( 9) and Eq. (27) in terms of the total mass M as 
seen by a distant observer 

IimM(r)=M. (32) 
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For the function tp (r) we shall always adopt the natural condition 



lim V (r) = 0. (33) 

r — >oo 

Inspection of Eqs. (25) and (27) reveals their different status. Indeed, Eq. (27) can 
easily be integrated for a general function M (r) and the final solutions is to be obtained by 
differentiation of the zeroth-order solution and making use of the boundary conditions. On 
the other hand, the first integral of the differential equation for M\ (r) cannot be constructed 
and one has to know the zeroth-order solution to determine M\. 

The assumed expansions of the functions M (r) and tp ( r ) as given by Eqs. (23) and (24), 
respectively, suggests that one can rewrite the boundary conditions of the first type in the 
following form: 

M (r + ) = ^±, Mi(r+) = 0, (oo) = 0, (34) 

whereas for the boundary conditions of the second type one has 

M {oo)=M, Mi(oo)=0, Vi(oo) = 0. (35) 

Now, let us concentrate on the zeroth-order equations supplemented with the conditions 
of the first type. Putting e = in Eq. (25), form (22) and (18) one obtains 

^ = ^fl-tanh^V (36) 
dr 2r 2 \ 26r / v ; 

which can be easily integrated to yied 

O 2 

M (r) = -6 tanh g- + C 2 . (37) 

Finally, making use of the conditions (34) one arrives at the desired result 

r_L O 2 O 2 

M (r) = ^- + b tanh b tanh ^- . (38) 

uw 2 2br + 2br y 1 

The thus obtained solution reduces to the Schwarzschild solution for Q = and it can be 
easily demonstrated that, by (26)and the boundary conditions (34) it remains so in the 
higher-order calculations. 

To specify the solution further we shall make use of the well-known relation [45] 

M = ^ - / s (2T b a TS b a ) K a dn, (39) 

where X is a constant time hypersurface and K a is a timelike Killing vector and apply it to 
the zeroth-order solution. Making use of the explicit form of the stress-energy tensor of the 
nonlinear electrodynamics one obtains 

M H = ^±+6tanh-^— , (40) 
2 2br + 

where Mh is the mass connected with the zeroth-order solution. We shall refer to Mh as to 
the horizon defined mass of the black hole. 

To develop the model further one has to determine the free parameter b. Our choice, 
which guarantees regularity of the zeroth-order line element at the center, is b = Mh, and 
hence Eq. (38) becomes 

Mo (r) =M H U- tanh . (41) 
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Unfortunately, the regularity of the zeroth-order solution docs not guarantee regularity of 
the higher-order perturbative solutions [30] . 

It should be noted that the Mr — Mh (Q, r+) and for fixed Q and r+ one has to determine 
Mr numerically. On the other hand it is possible to employ the equation M (r+) = r + /2 in 
the zeroth-order calculations and express all the results in (Q, Mh) parametrization instead 
of (Q,r + ). One can, therefore, construct solutions of this equation in terms of the Lambert 
function. Simple manipulations yield 



AM H Q 2 



AW+ {-peP) M H — Q 



2 ' 



where W + is a principal branch of the Lambert function and p 
solution for the inner horizon can be written in the form 

AM H Q 2 



AW- (-peP) M 2 H 



(42) 

Q 2 /4Mfj. Analogous 



(43) 



where W- is the second real branch of the Lambert function. (In fact, W + and W- are the 
only real branches.) Making use of the elementary properties of the Lambert functions one 
can demonstrate that the principal branch has the expansion 



W+{x) =x-x 2 + ^-x 3 - °x 4 + 0{x 5 ). 



(44) 



On the other hand, W-(x) — > — oo as x — > 0, and, consquently, the location of the event 
horizon tends to the Schwarzschild value whereas r_ — > 0. 

A typical run of Mjj as a function of £ for a few exemplary values of Q is showed in Fig. 1. 
For a given Q a line of Afg = const, intersects Q = const, curve at one or two points or it has 
no intersection points at all. The smaller one gives location of the inner horizon whereas the 
greater is to be identified with the event horizon. The minimum of Mh — Mh (Q = const, £) 
function represents extremal configurations when the two horizons merge. 




Figure 1: This graph shows solutions of the equation Mh = § + Mh tanh 2 m h g ^ or a ^ ew 
exemplary values of the charge Q. From bottom to top the curves correspond to Q = O.li, for 
i = 1, 12. For Mh = const., the greater solution represents location of the event horizon, 
r + whereas the smaller one represents the inner horizon, r_. The minimum of each curve 
corresponds to the extremal configuration with r + = r_ . 

It should be noted that the mass Mh is not the mass that would be measured for the 
perturbed black hole by an observer at infinity. Indeed, even for the zeroth-order solutions 
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the meaning of Mh and M is different and the substantial differences are transparent in the 
first order calculations. This can be easily seen by studying the limit 

M = lim M (r) = M H + eM x (oo) (45) 

r — >oo 

and Eq. (32). Identical result can be obtained form Eq. (39). Indeed, in order to apply (39) 
for the perturbed black hole one has to move the higher curvature terms of Eq. (9 ) into 
its right hand side and treat them as a contribution to the total stress-energy tensor. It can 
be demonstrated explicitly, that making use of Eq. (39) in the first-order calculations one 
obtains precisely (45). 

The function M\{r) can be expressed in terms of the polylogarithms. Unfortunately, 
it is rather complicated and to avoid unnecessary proliferation of long formulas it will not 
be displayed here. The first-order solution can be constructed employing the algorithm 
presented in Appendix of Ref. [30]. It should be noted that the function Mx(r) presented 
in [30] is calculated for the boundary conditions of the second type. 



4 The entropy 

Now, let us return to our main theme and calculate the entropy of the ABGB black hole. 
In doing so we shall put special emphasis on comparison of the results constructed for the 
nonlinear black hole with the analogous results obtained for the Reissner-Nordstrom solution. 
Such a comparison is especially interesting as the geometries of their classical counterparts 
are practically indistinguishable in two important regimes. To demonstrate this it suffices to 
expand the metric potentials in powers of \Q\/r+ and r+/r, respectively. Since the expansion 
takes the form 

2M H Q 2 Q 6 
f( r ) = 1 - — + -ZT - T77T72-J + ■■■> ( 46 ) 



12M 



2 r i 
H' 



the differences in the metric structure between ABGB and RN geometries in the exterior 
region for |Q|/r+ <C 1 arc small indeed. One has a similar behaviour for any (allowable) 
value of the charge for r»r + . 

The higher curvature terms in the action functional lead to the appearance of additional 
terms in the final expression describing entropy, which spoil area/entropy relation. Simple 
calculations carried out within the Noether charge framework indicate that the contribution 
of the quadratic part of the action to the entropy is given by 



SS = 2nr 2 + 



aR+±f3(Rl + R r r ) 



(47) 



Now, substituting the line element (13) with (14) and ( 23), into the general expression (47), 
expanding the right hand side of Eqs. (9) with respect to e, and, finally, retaining the linear 
terms only, one gets 



S = -nr\ + 2ner 2 



>^ + ) + ^M» ( , +) 
r + r + 



+ 0(s 2 ). 



(48) 



For the noncxtrcme black hole with the boundary conditions of the first type (34) one has 



2 2nQ 2 2 
6 =7rr + + — — —e cosh 



( Q 2 



M H r% 



\2M H r+ 



olQ tanh 



Q 2 



2M H r+ 



-(3 



Mnr+ — tanh 



( Q 2 



2 " " \2M H r+ 



)]} 



(49) 
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where Mh = Mh (Q, r+). In this approach the zeroth-order solution (41) determines the first 
order correction to the entropy of the nonextreme black hole completely. Having established 
Mh for given Q and r+ one can rewrite Eq. (49) putting S = S/M H , q = \Q\/Mh, 
x + = r + /Mu , a = a/M H and 8 = (3/M H . Simple manipulations yield 



TlX, + £ 



2irq 2 



cosh 



-2 Jt_ 

2x A 



aQ 2 tanh — 8 ( x + — — tanh — — 

2x+ V 2 2x+ 



(50) 



This results can be contrasted with the analogous result constructed for the Reissner-Nordstrom 
black hole 



c 2 »nQ 
S = nr + - 2(3 — — 



or 



S = tt(1 + V^f f - 2/3- 



irq 



(51) 



(52) 



To investigate the entropy S as given by Eq. (49) let us observe that for \Q\ /r + <C 1 one 
has r+ w 2Mh- Now, expanding hyperbolic functions in powers of \Q\/r + one obtains 



S = Trr z + - 2^8\ + O ( — ) 



(53) 



A comparison of Eqs. (53) and (51) shows that for \Q\ /r+ <C 1 the entropies of the ABGB and 
RN black holes are almost indistinguishable, as expected. It should be noted that contrary 
to the Reissner-Norstrom geometry, the entropy of the ABGB black hole depends on a and 
for |Q| /r + <C 1 the leading behaviour of this terms is ~ (Q/r + ) 4 . 

The analysis of the extremal configuration is more involved. First, let us return to the 
zeroth-order solution. It should be emphasized that although we do not ascribe any particular 
meaning to the zeroth-order solution, some of its features do possess clear and unambiguous 
meaning. For the boundary conditions (34) such a solution is described by the exact r + 
and Q. The extremality condition places additional relation between the elements of the pair 
(Q,r + ) or (Q,A4). Here we shall confine ourselves to the first pair. Simple considerations 
yield 

\Q\ = 2w 1 ' 2 M H (54) 



and 



Aw 
1 + w' 



M 



H, 



(55) 



where w = W + (l/e), and consequently 



\Q\ 



2W 1 / 2 ' 



(56) 



Returning to the first-order solution we recall the relation valid for the extremal configuration 
in the Reissner-Nordstrom geometry 

r+ = \Q\. (57) 

In Rcf. [12] we have ascribed this simple relation to tracelessness of the stress-energy tensor 
of the matter fields. As the stress-energy tensor of the nonlinear electrodynamics considered 
in this paper has a nonzero trace, one expects that the analogous relation between Q and r + 
in the ABGB geometry is more complicated. Indeed, after some algebra, one has 



2w 1 ' 2 \Q\ 
{1 + w) 



1 + e 



(3 + 2a 
16Q 2 w 



(w + 3) (w 2 - 1) 



(58) 
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Now, making use of (58) in (49) gives 



4 " " [{2a + f3)w 2 + 2(2a-f3)w + 5a + 2(3] , (59) 



eX " (l + w) 2 2{l + w) 

and the first term of the right hand side coincides with the Bekenstein- Hawking entropy [46] . 
Numerically, one has 

Sextr = kQ 2 x 0.6815 - 2tt£ x (0.0324a + 0.1754/3) , (60) 

where a common factor 2n has been singled out for convenience. Analogous relation for the 
extremal Reissner-Nordstrom black hole reads 

Sextr = kQ 2 - 27re/3. (61) 

Now, let us calculate the entropy of the ABGB black hole employing the Euclidean tech- 
niques propounded by Visser. First, observe that if the Lagrangian is arbitrary function of 
the Ricmann tensor (and its contractions) but is independent of its covariant derivatives, 
both methods, i. c. Wald's approach and Visser's method are equivalent. One may wonder, 
therefore, why we intend to carry out such a calculation. The answer is simple: although 
both methods should yield identical results, the calculational steps necessary to obtain the 
final result are quite different and consequently one can consider the calculations carried out 
within the framework of the one method as the useful check of the other. It is especially 
important in situations when the computational complexity of the considered problem may 
lead to numerous errors. 

The calculations proceed in a few steps. First, incorporate the Euclidean action functional 
of the quadratic gravity into the matter part of the action. Similarly, the (Lorentzian) energy 
density is given by 

P = -T} = T^£{F)-e(^-ll + i-Jl\. (62) 
H * 16tt v ' \8tt 1 8i V 

It could easily be demonstrated that pl — Le ~ O (s) and consequently it suffices to know 
the Hawking temperature to the zeroth-order. Moreover, due to subtle cancellations in the 
integrand of Eq. (6) the final result of the quadratures does not contain polylogarithm 
functions. Now, substituting 

i L Q 2 Q 2 



and (62) into (6), after some algebra, one has 



e 



5S = — t — r^K + N, (64) 



where rj = exp (Q 2 /2M H r+) 



4Q 6 20Q 4 \ 4 /20Q 4 72Q 4 12Q 



T -TK7 1*1 " -U7. — + ^7W + I *> + 



r + M§ M H ) \ M H r + r + M| rf M H 

(2QQ 4 12Q 6 56g 6 \ 2 ( 4Q 6 72Q 4 20Q 4 16g 6 , 

V + —772 - — + -77- - ZF77- V (65) 



\M H r+M* rfMu J \r+ M* r+ M H rfM H 
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and 



4Q 6 12Q 4 9 6Q 6 36Q 4 . , 

r^M H M H r + Mft r + ' 



V M H r+ M§ r+ r^M H J 

M H r+M% r+ rfM H J 
At first glance this result does not resemble Eq. (49). However, making use of the identity 

1-^-1, (67) 
r+ 

one can easily demonstrate that Eqs. (64-66) reduce precisely to Eq. (49). 



5 Final remarks 

In this paper we have constructed the entropy of the nonlinear ABGB-type black holes using 
the boundary conditions (34). The zeroth-order solution coincides, as expected, with the 
ABGB line element whereas the first-order correction can be elegantly expressed in terms 
of the polylogarithm functions. Now, let us briefly discuss the consequences of the second 
choice, in which the results are expressed in terms of the total mass of the system as measured 
by a distant observer. To calculate the location of the event horizon to the required order in 
(Q, M) parametrization one has to solve the first-order equations for Mi (r) and ipi (r) , and, 
subsequently, perturbatively solve the equation g tt (r + ) = assuming that the event horizon 
can be expanded as 

r+ =r a +er 1 +0(e 2 ). (68) 

Unfortunately, the function Mi (r) is rather complicated (it can be expressed in terms of the 
polylogarithms) and, once again, to avoid unnecessary proliferation of long formulas it will 
not be presented here. Interested reader is referred to [30]. 
Generally, for the nonexterme black hole one has 



S = 7T7-Q + 2nrie + 327tCTq 



—K (ro) + —r^K M 



+ O (e 2 ) 



(69) 



On the other hand, making use of (68), the equation (69) can be rewritten in the equivalent 
form 



S = nrl + e 



l-2M^(ro) +32 " CT ° 



4a , / x 2a + /?,,.. , . 
—Mq (r ) + —Mq (r ) 



- O (e 2 ) . (70) 



The extremal case should be analyzed separately. The extremal configuration of the 
ABGB black hole being the solution of the Einstein gravity is described by 



|Q C | = 2w 1 ' 2 M 



(71) 



and 



Aw 



M. 



(72) 
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One expects, that the higher-order curvature terms modify these relations, shifting (in a space 
of the parameters) extremal solution into a slightly different position. Indeed, treating Mq 
as a function of Q 2 and r, after some algebra, one concludes that the extremal configuration 
is still possible and is described by the relations 



Q 2 = Q 2 c +eA, r+=r c + s5, (73) 



where 



d_ 

dQ 2 



-l 



A = -(7^o) Ml (74) 



and 

6 = - (£ Mo ) 1 (l Mi ) + {w Mo & Ma ) 1 ( M w Mo ) • (75) 



Both S and A are to be calculated for the parameters describing extremal zeroth-order 
solution. Numerically, one has 

1.05314 0.43288^ 

and 

0.05121 0.57553 a 

s = -^T a -^^P- (77) 

Since the calculations of the entropy follow the general scheme sketched in previous section 
they will not be presented here. 

The purpose of the present paper (besides importance of the quadratic gravity in its own 
and the natural curiosity) is twofold. First, one can treat the calculations presented in this 
paper as the first step in understanding of the influence of the higher curvature terms on the 
entropy of black holes in a more complex setting than Maxwell electrodynamics. The next 
step would involve, for example, inclusion of the all curvature invariants of the order 4 and 6 
and degree 2 and 3 [47, 48, 49]. Moreover, it would be interesting to extend this analysis to 
general D-dimensional manifolds. The natural candidate for a higher-curvature theory is the 
Lovelock gravity [50]. Moreover, one may consider the more general curvature terms, with 
arbitrary coefficients rather than those inspired by particular theory. (See, for example [51, 47] 
and references cited therein.) On the other hand, and this is even more interesting, one can 
regard this sort of calculations as the preliminary results allowing to analyse and understand 
the typical subtleties one is likely to encounter when studying the semi-classical equations 
with the source term given by the renormalized stress-energy tensor of the quantized massive 
fields. Of course, the semi-classical equations are extremely complex [52, 8], but the general 
pattern that lies behind the calculations should remain the same. This group of problems 
are currently actively investigated and the results will be published elsewhere. 
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